Two-scale multitype contact process: coexistence 
in spatially explicit metapopulations 

Nicolas Lanchier 

Abstract It is known that the limiting behavior of the contact process strongly depends upon 
the geometry of the graph on which particles evolve: while the contact process on the regular 
lattice exhibits only two phases, the process on homogeneous trees exhibits an intermediate 
phase of weak survival. Similarly, we prove that the geometry of the graph can drastically 
affect the limiting behavior of multitype versions of the contact process. Namely, while it is 
strongly believed (and partly proved) that the coexistence region of the multitype contact 
process on the regular lattice reduces to a subset of the phase diagram with Lebesgue measure 
zero, we prove that the coexistence region of the process on a graph including two levels of 
interaction has a positive Lebesgue measure. The relevance of this multiscale spatial stochastic 
process as a model of disease dynamics is also discussed. 

1. Introduction 

The multitype contact process introduced in [15J is a continuous-time Markov process whose state 
space maps the d-dimensional integer lattice into the set {0,1,2} where state refers to empty 
sites and where state i, i = 1,2, refers to sites occupied by a type i particle. Denoting by r\t the 
state of the system at time t and by ~ the binary relation indicating that two vertices are nearest 
neighbors, the evolution of the process at vertex x is described by the transition rates 

c ^i(x,r/) = Bi ^ Hv(y) = !} cx-+o(x,rj) = 5 X 

x~y 

co^2{x,v) = B 2 ^2 1 i r l(,y) = 2 } c 2 -*o{x,rj) = S 2 

x~y 

where j(x, r/) is the rate at which the state of x flips from i to j. That is, type i particles give 
birth through the edges of the lattice to particles of their own type at rate Bi and die spontaneously 
at death rate <5j. If an offspring is sent to a site already occupied, the birth is suppressed. 

The multitype contact process has been introduced and completely studied when the death rates 
are equal by Neuhauser [15]. To fix the time scale, assume that 8\ = 8% = 1, and to leave out 
trivialities, assume in addition that the birth rates are greater than the critical value of the basic 
contact process |12| . Then, the type with the highest birth rate outcompetes the other type. In the 
neutral case when the birth rates are equal, the process clusters in dimension d < 2 while coexistence 
occurs in dimension d > 3. Here and after, coexistence means strong coexistence: there exists a 
stationary distribution with a positive density of type 1 and type 2. The long-term behavior of the 
process when the death rates are different remains an open problem but Neuhauser conjectured 
that her results extend to the general case provided one replaces the birth rate by the ratio of the 
birth rate to the death rate. In particular, it is believed that the coexistence region as a subset 
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of the space of the parameters has Lebesgue measure zero. The existence of such a coexistence 
region is mathematically interesting but for obvious reasons it is irrelevant to explain why species 
coexist in nature. In order to identify mechanisms (more meaningful for biologists) that promote 
coexistence, recent studies have focused on modifications of the multitype contact process in which 
the coexistence region contains an open set of the parameters. It has been proved in different 
contexts that coexistence is promoted by spatial [3J [J3] and temporal 0] heterogeneities. This 
article introduces the first example of a multitype contact process in which coexistence is produced 
by the geometry of the graph on which particles evolve. In some sense, our main result is analogous 
to the one of Pemantle [16] which states that, in contrast with the contact process on the regular 
lattice, the contact process on homogeneous trees exhibits a phase of weak survival. In both cases, 
the geometry of the graph is responsible for creating new qualitative behaviors. 

To construct our process, we consider the d-dimensional lattice Z rf as a homogeneous graph with 
degree 2d where vertices are connected to each of their 2c? nearest neighbors. Let N be an odd 
positive integer. Then, we consider the following collection of hyperplanes: 

H(i,j) = {x G M. d : x i = N/2 + jN} for i = 1,2, . . . , d and jeZ 

where x % denotes the ith coordinate of x, and remove from the original graph all the edges that 
intersect one of these hyperplanes. This induces a partition of the lattice into cf-dimensional cubes 
with length edge ./V that we call patches. See the left-hand side of Figure [U for an illustration where 
edges drawn in dotted lines are the edges to be removed. Since the parameter N is odd, each patch 
has a central vertex. To complete the construction, we draw a long edge between the centers of 
adjacent patches as indicated in the right-hand side of Figure [TJ The resulting graph can be seen 
as the superposition of two lattices that we call microscopic and mesoscopic lattices. Even though, 
for more convenience, we will prove all our results for this particular graph, our main coexistence 
result can be easily extended to more general graphs that we shall call two-scale graphs. These 
graphs are described in details at the end of this section. 

To formulate the evolution rules, we write x ~ y to indicate that vertices x and y are connected 
by a short edge, and x «->• y to indicate that both vertices are connected by a long edge. The 
evolution at vertex x is then given by the following transition rates: 



CQ^x(x,ri) = B x l{l(l/) = l} + Pi J2 1 My) = l} a^o(x,7i) = 5x 
c ^2(x,r]) = B 2 Hrj(y) = 2} + fa ^ l{ V (y) = 2} c 2 -+o(x,r}) = 5 2 . 



x^y xr^y 

Note that in the expression of Co-*. i{x,rf) the first sum is empty whenever vertex x is not located 
at the center of a patch. We call this process the two-scale multitype contact process. Also, we 
call the parameters (3t and Bi the microscopic and mesoscopic birth rates. The one-color version 
of this process has been introduced by Belhadji and Lanchier pQ as a spatially explicit model of 
metapopulation [101 H4|. The objective there was to determine parameter values for which survival 
occurs. In contrast, the emphasis here is on whether both types coexist or one type outcompetes 
the other type. While our analysis of the single-species model in [T] did not reveal any major 
difference between the contact processes on regular lattices and the graph of Figure HJ our analysis 
of the multispecies model shows that two-scale graphs, as opposed to regular lattices, promote the 
coexistence of the species. From now on, we assume that the parameters are chosen in such a way 
that each type survives in the absence of the other one and refer to [1] for explicit conditions of 
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Figure 1. Construction of the two-scale graph when N = 3 and d = 2. 

survival. Finally, note that when N = 1, i.e., patches reduce to a single vertex, the values of the 
microscopic birth rates are irrelevant and the two-scale multitype contact process reduces to the 
multitype contact process. Therefore, to avoid trivialities, we also assume that N ^ 1. 

We call the best invader the type with the highest Bi to 5i ratio, and the best competitor the 
type with the highest /3, to 5i ratio. Our first theorem extends Neuhauser's result to the two-scale 
multitype contact process: assuming that the death rates are equal, when one type is both the best 
invader and the best competitor, it outcompetes the other type except in the neutral case when 
the process clusters in d < 2 and coexistence occurs in d > 3. See Figure [2] for pictures of numerical 
simulations in the two dimensional neutral case. 



Theorem 1 (5 1 = 5 2 = 1 and JV / 1) Assume that B 1 < B 2 and fi\<fo- 
1. In the neutral case B\ = B 2 and f3\ = f3 2 we have the following alternative, 
(a) In d < 2 clustering occurs, i.e., for any initial configuration, 



lim P (jit(x) 

t— >oo 



1 and rjt{y) = 2) = for all x,y £ 



(b) In d > 3 coexistence occurs, i.e., there exists a stationary distribution under which the 
density of type 1 and the density of type 2 are both positive. 

2. If B\ < B 2 and (5\ < (5 2 with at least one strict inequality then type 2 wins. 

To search for strategies promoting coexistence, we now assume that one type, say type 1, is the 
best invader, and the other type is the best competitor, in which case the limiting behavior of the 
process is more difficult to predict. Interestingly, while the results of Theorem [1] are not sensitive to 
the patch size, the long-term behavior of the process under these new assumptions strongly depends 
upon the parameter N. As mentioned above, when patches reduce to a single vertex, the values 
of the microscopic birth rates are irrelevant so that type 1 particles outcompete type 2 particles 
as predicted by Theorem 1 in |15j . In contrast, taking iV large leaves enough room for type 2 to 
outcompete locally type 1 within each patch, except maybe near the central vertices located on 
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Figure 2. Snapshots of the multitype two-scale contact process on the 400 x 400 lattice with periodic boundary 
conditions at times 100 and 500, respectively. White and black vertices refer to species 1 and 2, respectively, and grey 
vertices to empty sites. The parameters are /9« = 3 and Bi = Si = 1 for i — 1,2. 

the mesoscopic lattice. The time a colony of type 2 particles persists within a single patch is long 
enough so that survival of type 2 is insured by casual migrations from one patch to another, which 
is referred to as the rescue effect in metapopulation theory. In conclusion, large patches promote 
survival of the best competitor, as indicated in the following theorem. 

Theorem 2 (5\ = 82 = 1) Assume that B\ > B2 > and fa > fii > (3 C . Then, in any dimension, 
type 2 survives provided the spatial scale N is sufficiently large. 

Theorem [5] is the key result to identify a set of parameters for which coexistence occurs. First 
of all, we fix the parameter values to make type 1 a good invader but a bad competitor living 
at a slow time scale. Type 1 then survives by jumping from patch to patch, or equivalently by 
invading the mesoscopic lattice. This holds regardless of the patch size. Coexistence of both types 
then follows from the proof of Theorem [2] by fixing the remaining parameters to make type 2 a 
good competitor living at a much faster time scale than type 1 and taking N large (which does 
not affect survival of type 2). To prove rigorously survival of both types, the two-scale multitype 
contact process will be simultaneously coupled with two different oriented percolation processes, 
one following the evolution of type 1 at a certain time scale, the other one following the evolution 
of type 2 at a slower time scale. The use of a block construction implies the existence of an open 
set of the parameters in which coexistence occurs, so we can conclude that 

Theorem 3 In any dimension, the Lebesgue measure of the coexistence region is strictly positive 
provided the spatial scale N is sufficiently large. 

As previously mentioned, Theorem [3] together with Neuhauser's conjecture |15j indicates that, in 
contrast to the regular lattice, the graph of Figure Q] promotes coexistence for the multitype contact 
process. Thinking of the competitive exclusion principle in ecology (the number of coexisting species 
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at equilibrium cannot exceed the number of resources), this suggests that the two-scale graph 
provides two spatial resources, namely the microscopic and the mesoscopic lattices, which allows 
two types to coexist. Note also that, while the coexistence region of the Neuhauser's competing 
model corresponds to the neutral case, for the two-scale multitype contact process, the coexistence 
region contains cases in which type 1 and type 2 have opposite strategies, namely one type is a 
good invader exploiting the mesoscopic lattice, while the other type is a good competitor using the 
microscopic lattice as its primary resource, and both types live at different time scales. 

Finally, even through for simplicity we will prove Theorems [2] and [3] only for the two-scale graph 
depicted in Figure[[J we would like to point out that our proofs easily extend to more general graphs, 
which also gives rise to realistic stochastic spatial models of disease dynamics. We first describe the 
general mathematical framework in which our results can be extended and then discuss about the 
relevance of this framework from a biological point of view. To begin with, let 

Hi = (Vi.Ei) and H 2 = (V 2 , E 2 ) with V x D V 2 and E x n E 2 = 

be two infinite graphs. We call H\ the microscopic graph and H 2 the mesoscopic graph, and consider 
the two-scale contact process evolving on the graph G = (Vi, E\ U E 2 ) where species i gives birth 
through the edges of the microscopic graph at rate and through the edges of the mesoscopic 
graph at rate -Bj. Assume that we have the following property that we call separation of the space 
scales. The mesoscopic graph H 2 contains a self-avoiding path {Xj : j G Z} C V 2 such that 

1. For all j £ Z, H x contains a self-avoiding path of length at least N containing Xj. 

2. For all x,y £ V 2 , the shortest path in H\ connecting x and y has length at least N. 

If there is no path in H\ connecting x and y (note that this is the case for the two-scale graph 
depicted in Figure [1]) we assume by convention that both vertices are connected by a path of infinite 
length so condition 2 above holds. Then, our proof of Theorem [3] implies that, for the two-scale 
multitype contact process evolving on G, the Lebesgue measure of the coexistence region is strictly 
positive provided N is sufficiently large. Note that, when H\ is a connected graph, condition 1 
above is always satisfied. In particular, a natural way to construct a suitable graph G is to start 
from an infinite connected graph Hi, then select an infinite subset V 2 C V\ of vertices that are at 
least distance N from each other, and finally add enough edges between the vertices in V 2 to obtain 
a mesoscopic graph H 2 with at least one infinite self-avoiding path. 

Returning to a microscopic graph made of infinitely many finite connected components and 
thinking of each component as a patch, one can legitimately argue that if two patches are connected 
by an edge of the mesoscopic structure then all the vertices of one patch should be connected to all 
the vertices of the other patch by a mesoscopic edge. In a number of contexts, however, patches are 
arbitrarily large yet adjacent patches are connected through only few vertices, in which case our 
general framework can capture the main features of the dynamics. This is the case, for instance, 
in metastatic diseases, such as malignant tumor cells that first spread within a given organ for a 
long time and then infect quickly other organs while reaching the bloodstream. In this context, 
the connected components of the microscopic structure represent organs or parts of the organs of 
the human body and the mesoscopic structure the vascular system. In a different context, one can 
think of the microscopic structure as a set of major cities and the mesoscopic structure as an airline 
network, where two types of diseases spread: one highly infectious disease such as H1N1 influenza 
that spreads quickly through the microscopic structure but slowly through the mesoscopic one due 
to airport screenings (best competitor), and one moderately infectious disease that spreads at an 
equal speed through the whole structure of the network (best invader). 
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rate 


symbols 


defined for each . . . 


effect on the configuration of the process 


Si 


x — > y 


edge (x,y) with x y 


if y is empty, it becomes of the same 
type as vertex x. 




2 

7* ^ 7/ 


edge (x, y) with x y 


if x is of type 2 and y is empty, y be- 
comes occupied by a type 2 particle. 


pi 


x- — >• J/ 


edge (a;, y) with x ~ ?/ 


if y is empty, it becomes of the same 
type as vertex x. 


pi -Pi 


2 

s — ► y 


edge (x, y) with x ~ y 


if x is of type 2 and y is empty, y be- 
comes occupied by a type 2 particle. 


1 


x at vertex x 


vertex x g 1 d 


if it exists, the particle at vertex x is 
killed regardless of its type. 



Table f 
Harris' graphical representation 



2. Proof of Theorem Q] 

Theorem [I] has been proved by Neuhauser [15] for the multitype contact process on the regular 
lattice, which corresponds to the case N = 1 for our process. Her proof relies on duality. Thinking 
of the process as being generated by a graphical representation, the dual process of the multitype 
contact process starting at a space-time point (x, T) exhibits a tree structure that induces an 
ancestor hierarchy in which the ancestors are arranged according to the order they determine the 
type of the particle at (x,T). Her result follows from the existence of a sequence of renewal points 
dividing the path of the first ancestor into independent and identically distributed pieces, as stated 
in details in Proposition 12.11 below. 

Similarly, the two-scale contact process is self-dual and the dual process exhibits a tree structure 
that allows us to identify a first ancestor. Renewal points can be defined from the topology of the 
dual process by using the same algorithm as for the multitype contact process. However, since the 
graph on which the particles compete is not homogeneous, the space-time displacements between 
consecutive renewal points are no longer identically distributed. The key to our proof is to rely 
on the fact that the graph of Figure Q] is invariant by translation of vector u G NIj d to show the 
existence of a subsequence of renewal points performing a random walk. 

To define the dual process, we first use an idea of Harris [TT] to construct the two-scale contact 
process graphically from collections of independent Poisson processes. These processes are defined 
for each directed edge (x, y) or vertex x as indicated in Table [H The last two columns show 
the rate of the Poisson processes and the symbols used to construct the graphical representation, 
respectively. Unlabeled arrows from x to y indicate birth events: provided site x is occupied and 
site y empty, y becomes occupied by a particle of the same type as the one at x. The same holds 
for type 2 arrows if the particle at site x is of type 2, but these arrows are forbidden for type 1 
particles, which takes into account the selective advantage of type 2. Finally, a x at site x indicates 
that a particle of either type at this site is killed. This graphical representation allows to construct 
the two-scale multitype contact process starting from any initial configuration. 

The construction of the dual process is done from the graphical representation by ignoring the 
labels on the arrows. We say that there is a path from (y, T — s) to (x, T), which corresponds to a 
dual path from (x, T) to (y, T — s), if there are sequences of times and sites 

so = T — s < 81 < — - < s n+ \ = T and xo = y,xi,...,x n = x 

such that the following two conditions hold: 
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(x,T) 

□ 



(2) 



(I) 



(2, 1, 1) 



(2, 2, 1, 2) 
(2,2, 1,2, I) 



(2, 2, 1, 1) 



Figure 3. Labeled tree structure of the dual process (d = 1 and N = 5). The bold lines refer to the path of the 
first ancestor which is determined by keeping track of the branch with the largest label. Numbers at the bottom of 
the structure indicate the ancestor hierarchy at a given time. White squares represent the sequence of renewal points 
assuming that the tree starting at the last white square at the bottom of the picture is infinite. 



1. For i = 1, 2, . . . , n, there is an arrow from Xj_i to Xj at time Sj and 

2. For i = 0, 1, . . . , n, the vertical segments {xj} x (sj, Sj+i) do not contain any x 's. 
The dual process starting at (x, T) is the set- valued process defined by 

f) s (x, T) = {y £ Z d : there is a dual path from (x, T) to (y, T — s)}. 

The dual process is naturally defined only for < s < T. Nevertheless, it is convenient to assume 
that the Poisson processes in Table Q] are defined for negative times so that the dual process is 
defined for all s > 0. The reason for introducing the dual process is that it allows us to deduce the 
presence of a particle at site x at time T from the configuration at earlier times by keeping track 
of potential ancestors. This appears in the duality relationship 

J]t(x) ^ <=^- VT-s(y) / for some y£r) s (x,T). 

The type of the particle at vertex x at time T can also be determined from the dual process and 
the configuration at earlier times: the dual process exhibits a tree structure that induces a so- 
called ancestor hierarchy. The members of the dual process at a fixed time (ancestors) are arranged 
according to the order they determine the type of (x,T). 

The ancestor hierarchy has been described algorithmically by Neuhauser. See the beginning of 
Section 2 in [15J. The original version of her algorithm is quite intuitive and, although it has been 
applied in a number of articles, did not evolve since then. We take advantage of this article to 



cS 



Nicolas Lanchier 



define rigorously the ancestor hierarchy from the tree structure of the dual process, following and 
improving an idea introduced in [T3]. Even if the dual processes of the contact process and two-scale 
contact process are different, their topologies are similar enough so that our approach applies to 
both models. The idea is to define a function <p s that maps the dual process at dual time s into 
the set of sequences with values in {1, 2, . . . } U {00} equipped with the lexicographic order <C- This 
function is strictly monotone in the sense that for xi,X2 £ fj s (x,T), 

4> s {xi) "C <As(x 2 ) X2 comes before x\ in the ancestor hierarchy. (1) 

The ancestor hierarchy on the dual process is thus naturally induced by the lexicographic order on 
the set of sequences through the function <f> s . Recall that 



(ui,U2,...) <C (vi,V2,...) if and only if 



Ui = Vi for i = 1, 2, . . . , n — 1 



for some positive integer n. The function cp s corresponds to a labeling of the tree structure of the 
dual process which is defined inductively by going backwards in time. The original branch starting 
at site x at dual time is labeled (00, 00, . . .). Now, assume that a birth event from x\ to X2 occurs 
at dual time sq to include site x\ to the dual process, and let 

4>s {x2) = (ui,n 2 , . . . ,n n ,oo, . . .) with m < 00 for i = 1, 2, . . . , n, (2) 

denote the label on the parent branch (the one at site X2 at dual time so). Then, set 

4> so (xi) = (ui,u 2 , . . . ,u n ,m, 00, . . .) 

if the new branch at site x\ is the mth one (going backwards in time) originated from the branch 
at site x 2 . Also, (J) s (xi) = <j) so (xi) until a death mark x is encountered at site x\ when the site is 
removed from the dual process. We refer the reader to Figure [3] for an example of realization of the 
dual process together with the label function along the branches of the tree. To avoid cumbersome 
notations, the sequences are identified to finite dimensional vectors: 

(ui,U 2 , ■ ■ ■ , 1%, OO, • • •) = (ui,U 2 , . . . ,U n ) 

where Ui < 00 for i = 1, 2, . . . , n. The label function is ideally suited to keep track of the topology 
of the dual process and has a number of interesting properties. For instance, if ([2]) holds then all 
the branches in the sub-tree starting at (x 2 , T — sq) have a label of the form 

(ui,u 2 , ■ ■ ■ ,u n ,u n+ i, . . . ,u n+k ,oo, . . .) for some u n+ i, . . . ,u n+k < 00. 

In particular, given the labels of two branches, it is straightforward to deduce the label of their 
most recent common ancestor (the root of the minimal sub-tree they both belong to) by looking 
at the first coordinates they have in common. The form of the sequence in ([2]) also indicates that 
the length of the dual path (number of arrows) from (x, T) to (x 2 , T — sq) is equal to n, which can 
be seen as a number of generations. Also, it is easy to check that the description of the ancestor 
hierarchy given in Section 2 of [15] reduces to ([T|) above. 

The path of the first ancestor, which is the object of primary interest to determine the type of 
the particle at (x,T), is constructed by following backwards in time the branch with the largest 
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label as defined in condition (pQ). The path of the first ancestor is a complicated object. However, in 
the case of the multitype contact process, Neuhauser [JS] proved that this path can be divided into 
independent and identically distributed pieces at some particular space-time points called renewal 
points. To define renewal points, we say that (X, r) lives forever if 

t) s (X,t) ± for all s > 0. (3) 

Then, whenever the path of the first ancestor jumps to a point (X, r) that lives forever, this point 
is a renewal point. Note that renewal points are well defined only if the starting point (x, T) lives 
forever, which we assume from now on. Let {(X n ,T n ) : n > 0} denote the sequence of renewal 
points starting from (Xq, tq) = (x, 0). See Figure [3]for a picture. For the multitype contact process, 
which is the case N = 1 for our process, Neuhauser [15] proved the following proposition. 

Proposition 2.1 (Neuhauser, 1992) Assume that N = 1. Then {(X n ,T n ) : n > 0} perforins a 
random walk on Z d x M . . Moreover, there exist C\ < oo and 71 > such that 

P(\\X n+ i - X n \ > s) < Ci exp(-7is) and P (r n+1 - r n > s) < C\ exp(~7is) 

for all n > and all s > 0. 

Proposition 12.11 is the key to understand the ancestry of a single site and deduce how different 
sites are correlated. Unfortunately, it does not hold when N > 3 in which case the space-time 
displacements between consecutive renewal points are no longer identically distributed due to the 
geometry of the graph in Figure [TJ Observing that the parent of a particle located at the corner of a 
patch has to be closer to the center than the particle itself, we see that the path of the first ancestor 
undergoes a drift directed to the center of the patch. The drift is stronger while approaching the 
boundary of the patch. Also, the path of the first ancestor is more likely to jump from one patch to 
another while approaching the center of a patch. The key to prove Theorem Q] is to observe that the 
graph is invariant by translation of vector u G N1j d . This implies that the distributions of the dual 
processes starting at x and x + u, respectively, can be deduced from one another by a translation 
of vector ±u, and suggests the existence of subsequences of renewal points satisfying Proposition 
12.11 More precisely, we extract a subsequence inductively by letting 

no = and n^i = min{n > Hi : X n £ N7* d } for i = 0, 1,2, . . . 

Note that, regardless of X n , the probability that X n+ i G NIj d is bounded from below by a positive 
constant that only depends on N. This, together with the Borel-Cantelli Lemma and the fact 
that what happens before and after a certain renewal point uses disjoint parts of the graphical 
representation and so is independent, implies that ni is almost surely finite for all i. Let 

Yi = X ni and Oi = r ni for all i > 0. 

Then, {(Y^, crj) : i > 1} represents the subsequence of renewal points visiting the center of patches 
and we have the following proposition. 

Proposition 2.2 The subsequence {(Y{,ai) : i > 1} is a random walk on N7L d x M + regardless of 
the value of N > 1. Moreover, there exist C2 < 00 and 72 > such that 

P(\\Y i+ i - Yi\\ > s) < C 2 exp(-7 2 s) and P (a i+ i - Oi> s) < C 2 exp(-7 2 s) 

for all i > 1 and all s > 0. 
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Proof. Proposition 12.11 indicates that the increments (X n+ i,T n+ i) — (X n ,T n ) are independent for 
different values of n. This follows from the fact that, after a renewal point, the path of the first 
ancestor is determined by the structure of the sub-tree starting at the renewal point (rather than 
the "super structure" of the entire dual process) and so depends only on parts of the graphical 
representation that are after the renewal point. The same holds for {(li,<7j) : i > 0} because it is a 
subsequence of {(X n ,r n ) : n > 0}. Now, since the graphical representation is translation invariant 
in time and the graph of Figure [T] is invariant by translation of vector u G Nl, d , we have that for 
all s > and any collection A of subsets of 7L d 

P (fj a (x,T) - x G A) = P (f) a (y,S) - y G A) whenever x-yeNZ d . 

In particular, using that Y i+1 -Yi£ NZ d for i > 1, we obtain 

P(f) s (Yi,T- <n) -Yi€A) = P(f) a (Y i+1 ,T - a i+1 ) - Y i+1 G A) for all i > 1. 

Note however that this does not hold for % = since Yq = x is not a priori at the center of a patch. 
Since (Yj+i, fj+i) is determined from {fj s (Yi,T — ai) : s > 0}, we deduce that 

P{(Y i+ll a i+1 )-{Y h cn)eB) = P((Y i+2 ,a i+2 )-(Y i+1 ,a i+1 ) e B) 

for any measurable set B C Z, d x M + and all % > 1. In conclusion, the space-time displacements 
are identically distributed which makes {(1^, crj) : % > 1} a random walk. To prove the exponential 
bounds on the space-time displacements, we first observe that, the parameter ./V being fixed, there 
exists a positive constant > such that 

P (X n+1 G NZ d \X n = y) > p N for all y G Z d and all n > 0. 

In words, the next renewal point visiting the center of a patch can be found after at most a geometric 
number of steps, which can be expressed formally as follows: 

P (n^+i — rii > j) < P (K > j) for all j > 1 where K ~ Geometric (pn)- 

Now, members of {t„ + i — r n : n > 0} are clearly independent. Moreover, the same arguments as in 
the proof of Proposition 12.11 that can be found in [15] imply the existence of a continuous random 
variable r with exponentially bounded tails such that 

P (r n +i - r n > s) < P (t > s) for all s > and n > 0. 

Let M(X, 0) stand for the moment generating function of X at point 9. Since K is geometri- 
cally distributed and temporal increments between consecutive renewal points are independent and 
uniformly stochastically smaller than r, there exists 6q > 1 such that 

M(a i+1 -o-i,0) = M(T n . +1 -T m ,9) < M(t k ,0) 

oo 

< P(K = n) M(r 1 ,e)M(T 2 ,e) ••• M(r n ,9) 

n=l oo 

< P(K = n) [M(T,6)] n = M(K,M(t,9)) < oo 

71=1 
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for all 6 £ (1,6q). By using Markov's inequality, we obtain that for 9 £ (l,#o) fixed, 

P(a i+1 -Oi> s) < 6~ s M(a i+1 - a h 9) = C 2 exp(-7 2 s) 

where C 2 = M{ai + i — < oo and 72 = log# > 0. The exponential bound on the spatial 

displacements between consecutive renewal points follows from the exponential bound on the tem- 
poral displacements since the dual process grows at most linearly. This follows from the fact that 
the two-scale contact process is self-dual and that the lengths of all the invasion paths at a given 
time t are uniformly bounded stochastically by iV times the radius of the Richardson's model with 
parameter max(l?2,/32) at time t. □ 

With Proposition 12.21 in hands, Theorem [1] in the general case when N > 1 follows from the 
techniques developed in Sections 2-5 of [15]. The main idea in the neutral case is that the random 
walk {Yi : i > 1} is recurrent in d < 2 which implies that the paths of the first ancestors of a finite 
number of sites eventually coalesce with probability 1. This makes these sites identical by descent 
which translates into a clustering of the system. In contrast, transience of the random walks in 
dimension d > 3 implies that two distinct sites may not be identical by descent in which case the 
states of both sites are independent whenever all the sites are independent at time 0. 

3. The multitype contact process in finite volume. 

In preparation for proving Theorems [2] and O we investigate the restriction of the multitype contact 
process to a single patch, i.e., particles landed outside the patch are killed. The process is further 
modified to include spontaneous births of type 1 particles at the center of the patch. The reason for 
introducing spontaneous births is to later compare this process with the two-scale contact process 
viewed on a single patch in which offspring may originate from other patches. More precisely, the 
state of the process at time t is a spatial configuration 77^ : Aq — > {0, 1, 2} where 

A = (-N/2,N/2) d n Z d 

is the patch centered at 0. The evolution at x £ Aq is described by 

co^ 2 (x,v N ) = h Hv tf (y) = 2} 

ci^o(x,f] N ) = c 2 ->o(x,r) N ) = 1. 

where x ~ y indicates that \\x — y\ = 1 and x,y £ Aq. In order to understand both the multitype 
contact process in finite volume and the two-scale contact process, the first step is to compare the 
processes viewed on suitable length and time scales with oriented percolation. To begin with, we 
introduce 1-dependent oriented percolation with parameter 1 — e on 

g = {( z , n) £ Z d x Z + : z 1 + . . . + z d + n is even} 

where z % denotes the ith coordinate of z £ Z d . Each site (z, n) £ Q is associated with a Bernoulli 
random variable u(z,n) £ {0, 1} with success probability 1 — e, and is said to be closed if there is 
a failure (0) and open if there is a success (1). The 1-dependency means that 

P (co(zi,ni) = 1 for 1 < i < m) = (1 — e) m whenever \\(zi,rii) — (z,j,rij)\ > 1 for i ^ j. 
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A site (z, n) is said to be wet (at level n) if there exist zq, zi, . . . , z n = z such that 

1. For i = 0, 1, . . . , n — 1, we have f^+i — z,- L \ = 1. 

2. For i = 0, 1, . . . , n, site (zi, i) is open. 
From now on, the parameter e > is fixed so that 

P (W n + for all n > | W = {0}) > 0, 

where W n denotes the set of wet sites at level n. To investigate the interacting particle system 
restricted to a single patch, we let ./V = (K + 2)L with K and L odd, which induces a partition of 
patch Aq into L-cubes, and also consider oriented percolation on 

Gk = Q n {(-K/2,K/2) d x Z+} 

= {(z, n) G T, d x Z + : z 1 + . . . + z d + n is even and supj \z l \ < (K — l)/2}. 

For the oriented percolation process on Qk we denote by co K (z,n) the Bernoulli random variable 
associated to site (z, n) and by Wjf- the corresponding set of wet sites at level n. The proof of 
Theorems [2] and [3] relies on multiscale arguments. More precisely, we will consider the following 
three mesoscopic spatial scales: 

1. Upper mesoscopic scale. The lattice Z d is partitioned into iV-cubes called patches. In the next 
section, we compare the two-scale contact process viewed at the patch level with the oriented 
percolation process on Q introduced above. 

2. Intermediate mesoscopic scale. Each iV-cube is partitioned into L-cubes: 

B z = Lz + B where B = (-L/2, L/2) d n 7L d and z G Z d . 

In this section, we compare the multitype contact process on patch viewed at the L-cube 
level with the oriented percolation process on Qk introduced above. 

3. Lower mesoscopic scale. Each L-cube B z is further divided into L 01 -cubes: 

D w = L 0A w + D where D = (— L°' 1 /2, L°' 1 /2) d n Z d and w G 7L d . 

Since Aq is finite, type 2 particles go extinct eventually, and the process r)f converges weakly to 
a stationary distribution with a positive density of type 1. The main objective of this section is to 
prove that when /?2 > P\ and starting with a significant number of type 2 particles the time to 
extinction grows exponentially with the size of the patch Aq. Moreover, locally in space and time, 
type 2 outcompetes type 1 in the sense that, excluding a neighborhood of site whose size does not 
depend on N, most of the patch is void of l's and has a positive density of 2's up to the extinction 
time. To prove these results, the first step is to couple the process with the oriented percolation 
process on the lattice Qk through the following 

Definition 3.1 Let = (-L/6,L/6) d and T = I?. We call (z,n) G Q K a good site if 

1 . The set B z \ B* is void of 1 's at time nT and 

2. For all w such that D w C B Z \B*, the set D w contains at least one 2 at time nT. 
The set of good sites at level n is denoted by 



X% = {z G 1 d : (z, n) G Qk and (z, n) is good} 
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H pfi n pH frtT p a p n 






a — ► y 


edge (a;, y) £ Aq with a; ~ y 


if y is empty, it becomes of the 
same type as vertex x. 




2 

x — ► y 


edge (a;, y) 6 v4q with x ~ y 


if x is of type 2 and y is empty, y 
becomes occupied by a type 2 par- 
ticle. 


1 


x at vertex x 


vertex a; £ Z d 


if it exists, the particle at vertex x 
is killed regardless of its type. 


2dB! 


• at vertex 


vertex only 


if x is empty, it becomes occupied 
by a particle of type 1. 



Table 2 
Harris' graphical representation 



Note that site (0, n) has a special treatment: there is no requirement about the spatial configuration 
of the system inside B* C Bq. The reason is that, due to the presence of spontaneous births, there 
is a region around the center of the patch which is occupied by l's most of the time. Note also 
that we have excluded boxes B z with supj \z l \ = {K + l)/2 which are located along the frontier of 
patch Aq since, due to boundary effects, the density of type 2 particles in these boxes can shrink 
significantly. The key result of this section is the following 

Proposition 3.2 Let 02 > Pi > Pc- Then, for L = L(e) sufficiently large, the processes can be 
constructed on the same probability space in such a way that 

P (W£ C X* for all n > | W* = X?) = 1. 

In order to establish Proposition l3.2[ the objective is to prove that there exists L = L(e) sufficiently 
large such that the following holds: 

P ((z2,n + 1) is good | {z\, n) is good) > 1 — e 

for all z\, Z2 with \z\ — Z2I = 1 an d (zi,n), {z%^ n + 1) E Qk- 

With inequalities (j4]) in hands, the proof of Proposition 13.21 follows from standard techniques, and 
we refer the reader to the Appendix of [6] for more details. Since the central site plays a particular 
role due to the spontaneous births of type 1 particles as pointed out in Definition 13.14 there are three 
cases to be considered, as illustrated in Figure H] (the three pictures in this figure will be explained 
later). We only prove ([!]) in the case when Z2 = (left picture) which is slightly more complicated 
than the two other cases. Using spatial symmetry and homogeneity in time of the evolution rules 
of the process, it suffices to prove the following 

Lemma 3.3 Let e\ = (1, 0, . . . , 0) G 7L d . Then 

P ((0, 1) is good I (ei, 0) is good) > 1 — e for L sufficiently large. 

The proof of Lemma 13.31 relies on duality techniques as well. In order to define the dual process, 
the first step is to construct the multitype contact process on Aq graphically from collections of 
independent Poisson processes [11] . As indicated in Table [2j these processes are defined for each 
directed edge (x,y) G Aq x Aq or vertex x G Aq. Unlabeled arrows, type 2 arrows, and x's have the 
same interpretation as in Table [1] above. The additional symbol • indicates a spontaneous birth of 
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type 1 particle at site 0. This graphical representation allows us to construct the multitype contact 
process on Aq starting from any initial configuration. 

We say that there is a path from (y,T — s) to (x,T), or equivalently that there is a dual path 
from (x, T) to (y, T — s), if there are sequences of times and vertices 

so = T — s < si < ••• < s n+ i = T and xo = y,x\,...,x n = x 

such that the following two conditions hold: 

1. For i = 1, 2, . . . , n, there is an arrow from Xi-i to x% at time S{ and 

2. For i = 0, 1, . . . , n, the vertical segments {xi} x (sj, Si+i) do not contain any x 's. 

Note that, in our definition of path and dual path, «'s have no effect, though they are important in 
the construction of the process. The dual process starting at space-time point (x, T) is then defined 
as the set-valued process 

fj s (x, T) = {y G Z d : there is a dual path from (x, T) to (y, T — s)}. 

As previously, it is convenient to assume that the Poisson processes in the graphical representation 
are defined for negative times so that the dual process is defined for all s > 0. To deduce the 
type of the particle at (x, T) from the configuration at earlier times, we define a labeling of the 
tree structure of the dual process, thus inducing an ancestor hierarchy, by using the algorithm 
introduced in Section [2j The path of the first ancestor is then constructed by following backwards 
in time the branch with the largest label. The type of (x, T) is determined as follows: 

1. If the first ancestor crosses at least one • on its way up to (x, T) ... 

(a) regardless of the initial configuration, (x,T) is of type 1. 

2. If the first ancestor does not cross any »'s on its way up to (x, T) . . . 

(a) and lands at time on an empty site, the first ancestor does not determine (x,T). 

(b) and lands at time on a 1 and that, on its way up to (x, T), the first ancestor does not 
cross any type 2 arrow then (x, T) is of type 1. 

(c) and lands at time on a 1 and that, on its way up to (x,T), the first ancestor crosses 
a type 2 arrow then the first ancestor does not determine the type of (x,T). We then 
follow the path of the first ancestor on its way up to (x, T) until the first 2-arrow we 
encounter and discard all the ancestors of the point where this arrow is directed to. 

(d) and lands at time on a 2 then (x, T) is of type 2. 

If the first ancestor does not determine the type of (x,T) (2a and 2c above), we look at the next 
ancestor in the hierarchy, and so on. Point la follows from the fact that a spontaneous birth of type 
1 particle occurs along the path of the first ancestor. Points 2a-2d are the same as for the multitype 
contact process and we refer the reader to [15], page 472, for more details on how to determine the 
type of (x, T) for the process with no spontaneous birth. Finally, note that, since the state space 
of the process is finite, 

lim fj (X, r) = with probability 1 

which contradicts the definition of "living forever" introduced in ([3]) above. In this section, we say 
that a space-time point (X, r) lives forever if 

f?f(X,r) / for all s < r, (5) 
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that is, condition ([3]) is only satisfied for all s < r. Then, whenever the path of the first ancestor 
jumps to a space-time point that lives forever in the sense of ([S]), this point is called a renewal 
point. As in Section [2l the sequence of renewal points divides the path of the first ancestor into 
independent and identically distributed pieces. We are now ready to prove Lemma 13.31 

Proof of Lemma 13 . 31 Let x E Bo \ B* and assume that (x,T) lives forever. The first step is to 
prove that there exist C3 < 00 and 73 > such that, for all L sufficiently large, 

P(rj T (x) / 2 I (ei,0) is good and (x,T) lives forever) < C3 exp(— 73L ' 1 ). (6) 

To prove ([6]), we will construct a dual path & s forbidden for the l's starting at (x,T) and ending at 
time on a site occupied by a type 2 particle. The idea is to apply a modification of the so-called 
repositioning algorithm (for the original version, see page 28 in [8]). We say that a renewal point 
is associated with a 2-arrow if the first arrow a particle crosses starting at this renewal point and 
moving up the graphical representation is a 2-arrow. We call selected path 2l s with origin x and 
target y E Bq the following dual path: The process starts at 2lo = x and follows the path of the 
first ancestor starting at (x,T) until the first time a\ it jumps to a renewal point associated with 
a 2-arrow. Then, we either leave 2l s where it is at that time or reposition it. To determine whether 
and where to reposition the selected path, we denote the location of the second ancestor in the 
hierarchy at time a\, provided this ancestor exists, by 5S (J1 . Let m be a large constant that does not 
depend on L, and let A be the straight line going through x and y. Also, we denote the Euclidean 
distance by dist( • , • )• 

1. Assume that *8 CT1 exists and (tB ai ,T — o\) lives forever. Then 

(a) if dist(2l CTl , A) > m and dist(*B CTl , A) < dist(2l CTl , A) set 2l CTl+ = 5S CTl . 

(b) if dist(2l CT1 , A) < m and dist(*8 CT1 , y) < dist(2l CT1 , y) set 2l CT1+ = <B CT1 . 

(c) otherwise, we set 2l<j 1+ = 2l CTl . 

2. Assume that *B (Jl does not exist or (Q5 CTl ,T — a%) does not live forever. Then 

(a) we set 2lo-i+ = 21<ti • 

In either case, we start a new dual process at (2t CTl +,T — <t\) and follow the path of its first ancestor 
until the first time 02 it jumps to a renewal point associated with a 2-arrow when we apply again 
the repositioning algorithm, and so on. Intuitively, this causes the selected path 2l s to drift towards 
the target y while staying close to the straight line A. More precisely, let x\ and xi belong to the 
segment (x, y) in the order indicated in the following picture: 



Xl 



Xz 



A 
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Assume that \x\ — X2I = VL and ||2l So — xi|| < at some time so and set 

S Xl = inf {s > sq : dist(2l s ,xi) > VT} 
T X2 = inf {s > s : dist(a s ,x 2 ) < VT/4}. 

Then, it can be proved that, for suitable constants C^C*, < 00 and 75 > 0, 

P(S X1 < T X2 or T X2 > C A Vl) < C 5 exp(- 75 L a2 ). (7) 

The proof of j7|) can be found in [13], Lemma 3.5. That is, starting from the small ball centered 
at x\ in the picture above, the selected path hits the small ball centered at X2 before leaving the 
large ball centered at x\, this takes less than C^yL units of time with probability close to 1 when 
the parameter L is large. In particular, if we let 

S A = inf {s > : dist(2l s , A) > \[L or dist(2l s , 0) > 2L} 
T y = inf {s > : dist(2l s , y) < VZ/4} 

apply (|7|) consecutively and use that dist(x,y) < Ly/d, we obtain 

P (S A < T y or T y > C^LVd) < C^yfdyfl exp(-7 5 L a2 ). (8) 

To construct & s , we let y be the corner of B* = (—L/3,L/3) d closest to x and z = Le\ be the 
center of patch B ei (see Figured]). The dual path & s starts at ©o = x, follows the selected path 
with target y until T y when it hits the Euclidean ball with center y and radius y/L/4, then follows 
the selected path with target z until time r = T — \[L. By (jHJ) , 

P (dist(6 s , 0) < VT for some s < r or dist(6 T , z) > v 7 !) < C 6 exp(-7 6 L ' 2 ) (9) 

for suitable Cq < 00 and 7 6 > 0. In other respects, since (ei, 0) is good, a straightforward application 
of Lemma 3.9 in [13J implies that 

P ((@T) VL) is not occupied by a 2 | dist((5 T , z) < vL 

and (ei, 0) is good) < C7 exp(— 7 7L 01 ) 

for suitable C7 < 00 and 77 > 0. Finally, using the duality properties described above and the fact 
that 2-arrows are forbidden for the l's, and combining Q and (|l(jp . we obtain 

P(r)!p(x) ^ 2 I (ei,0) is good and (x,T) lives forever) 

< P ((& T , VL) is not occupied by a 2 | (ei, 0) is good) 

+ P(dist(6 s ,0) < VI for some s < r) < C 3 exp(-7 3 L 01 ) 

for suitable C3 < 00 and 73 > and all L sufficiently large, which establishes ([6]). The second step 
is to prove that 

P{rj^{x) = 1 I (x, T) does not live forever) < C$ exp(— 7§L) (11) 



for appropriate Cs < 00 and 78 > 0. Note that, on the event that (x,T) does not live forever, 
duality implies that the probability in (|lip is equal to for the multitype contact process with no 
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spontaneous birth of type 1 particles. In our case, due to the presence of spontaneous births at the 
central vertex 0, we need to bound the probability that the dual process starting at (x, T) hits 0. 
More precisely, since \\x — Of > L/6, (jlip follows from 

P (0 E f}^ (x, T) for some s < T \ rj^ (x, T) = 0) 

< P (the dual process starting at (x,T) has radius 

at least L/6 | (x,T) does not live forever) < Cg exp(— 7§L) 

which is a well-known property of the contact process. The proof follows from the analogous result 
for oriented percolation (see [5], Section 12) and the fact that the supercritical contact process 
viewed on suitable length and time scales dominates oriented percolation (see [2]). Finally, combin- 
ing © and (fTT|) . we obtain 

P (t]t (x) = 1 for some x E Bq\ | (ei, 0) is good) 

< C 3 L d exp(- 73 L - 1 ) + C 8 L d exp(- 78 L) < C 9 exp(- 79 L ai ) iU) 

for suitable Cg < oo and 79 > 0, and all L sufficiently large. That is, condition 1 in Definition 13.11 
holds with probability arbitrarily close to 1. Now, let D w C Bq \ B*. Since the process dominates 
a one-color contact process with parameter /3i > /3 C , we have 

P (t)t (x) = for all x E D w | (ei,0) is good) < C w exp(-7i L ai ) (13) 

for appropriate C\$ < 00 and 710 > 0. Combining (fl2|) and (fl3|) implies that 

p (r]N (3;) 2 for all x E D w \ (ei, 0) is good) 

< P (t]t (x) = 1 for some x E D w \ (ei,0) is good) 

+ P (r^ (x) = for all x E D w \ (e u 0) is good) 

< C 9 exp(- 79 L - 1 ) + C10 exp(- 7l0 L ai ). 

In particular, there exist C\\ < 00 and 711 > such that 

P (there exists D w C Bq \ B* : r/^(x) 7^ 2 for all x E D w \ (ei, 0) is good) 

< C 9 L d exp(- 79 L°- 1 ) + C7 10 L d exp(-7i L ' 1 ) < C u exp(- 7 nL - 1 ) U Jj 

for L sufficiently large. The lemma then follows from (|12p and ()14|) . □ 

From Lemma [3.3l it is easy to deduce (HJ) when z<i = 0. The left picture of Figure U] gives a schematic 
illustration of the dual path & s in continuous line. The proof of when z\ = and when both 
z\ and Z2 are 7^ is similar and we refer to the last two pictures in Figured] for an illustration of a 
suitable dual path in these two cases. Note that the case Z2 = is slightly more complicated since 
the repositioning algorithm has to be applied in two different directions in order to avoid with 
high probability. With Proposition 13.21 in hands, we are now ready to deduce useful properties of 
the multitype contact process restricted to Aq from analogous properties of the oriented percolation 
process on Qk- We let 

tat = inf {t > : rj? (x) / 2 for all x E All- 
denote the extinction time of type 2 particles. We denote by Pq the conditional probability given 
the event that vertex (and only vertex 0) is occupied by a type 2 particle at time 0. In the next 
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Figure 4. Picture of the selected path. 



lemma, we prove that, starting with a single 2 at the center of the patch, with high probability, 
the process may exhibit only two extreme behaviors: either the 2's spread out successfully and the 
time to extinction is arbitrarily large when ./V is large, or they die out quickly. 

Lemma 3.4 Let /3 2 > Pi > Pc an d %K = exp(c.K'). Then, there exists c > such that 

P (t <t n < 2>Z K ) < exp(-7i 2 K) + Ci 3 exp(-7i 3 t) 
for all K sufficiently large and suitable constants C12, C13 < 00 and 712, 713 > 0. 

Proof. The proof is divided into two steps, both relying on Proposition 13.21 The idea is to decom- 
pose the event of interest according to whether the event 

{Ik = {(Xn , n ) <t- Qk-2 for some n > 0} 

occurs or not. The occurrence of {Ik can be thought of as a successful invasion of type 2. We will 
prove that (i) type 2 particles live an exponentially long time on the event {Ik, while (ii) they die 
out quickly on the complement of {Ik- 

(i) The event {Ik occurs. In this case, we will prove that 

P ({Ik and t n < 31k) < Cn exp(-7i 2 K) (15) 

for suitable constants Cxi < °o and 712 > 0. Let 

n = min{n>0: (X* ,n) G K -2}- 

Since the event {Ik occurs and the set dominates the set of wet sites of a supercritical perco- 
lation process, there is an in-all-direction expanding region centered at which contains a positive 
density of good sites. Since the correlation between two sites decays exponentially with the distance, 
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the configurations in two Lx L squares with no corner in common are almost independent when the 
parameter L is large, which, together with large deviation estimates for the Binomial distribution, 
implies the existence of a constant 714 > such that, for K large, 

P (Q K and cardX,f < lu K) < C 15 exp(-j 15 K) (16) 

for suitable C15 < 00 and 715 > 0. Now, assume that card W^ > 714-^ and let 

tt k = inf{ra > : = 0} 

denote the extinction level of the percolation process restricted to Qk- Theorem 2 in [9] implies 
that there exists c > 0, fixed from now on, such that, for K sufficiently large, 

P(tt k < ST^exp^K) and cardJU^ > <y u K) < C m exp(- lm K) (17) 

for suitable Ci6 < 00 and 716 > 0. Combining estimates (fT6j) and (fT?]) with the coupling provided 
in Proposition 13.21 implies that 

Pq (Qk arid 77/ < 3Z^-) < Pq (Qk and cardX^ < J14K) 

+ P (t n < KL K and card X* > < Pq (Q k and card X* < juK) 

+ P(tt k < ZT- l l K and cardW^ > lu K) < C 15 exp(- 7l5 ^) + C 16 exp(- 7l6 K). 

This completes the proof of (fl5l) . 

(ii) The event Qk does not occur. In this case, we will prove that 

P (Q C K and r N >t) < C 13 exp(- 7l3 t) (18) 

for suitable constants C13 < 00 and 713 > 0. By Proposition 13.21 it suffices to prove the analogous 
result for the oriented percolation process, namely 

P ((W*, n) C Gk-2 for all n > and tt k > m \ W* = {0}) < C 17 exp(-j 17 m) (19) 

for some Cyj < 00 and 717 > where ttk is the extinction level. To establish (fl9|) . we couple the 
restricted and unrestricted percolation processes by letting 

W = W K = {0} and w(z,n) = co K (z,n) for all (z,n) G Qk- (20) 

Assume that there exists x G W n \ . This together with (|20p implies that any open path ending 
at site (x,n) for the unrestricted percolation process leaves the set Qk- Moving up along such a 
path, we denote by (xq^tiq) the first site outside Qk-2 we encounter. Using ([20]) again, it is easy 
to see that there is an open path ending at (xQ,no) for the restricted percolation process which 
implies that (W^,n) (jL Gk-2 for some n > 0. In conclusion, 

W n / for some n > => (W*,n) <£ Q K -2 for some n > 0. (21) 

Using the reverse of (|21|) . we can bound the left-hand side of (|19|) by 

P ((W n , n) C Qk-2 for all n > and vr > m \ W = {0}) 

< P (W n = for some n > and 7r > m \ Wo = {0}) 

= P (m < 7r < 00 I Wq = {0}) < C17 exp(— 717m) 
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where n = inf {n > : W n = 0} which is a well-known property of supercritical oriented percolation 
processes (see page 1031 in [5]). This establishes (fTUj) and fTS]) . 
We conclude by noticing that 

Po (t<T N < 3I K ) < P (Sl K and t n < 3I K ) + P (n c K and t n > t) 
< Ci2 exp(-7i2-K') + C13 exp(-7i 3 t). 

This completes the proof. □ 

Due to spontaneous births, the region near vertex is mostly occupied by l's. The next two lemmas 
show however that, provided the 2's invade the patch successfully, the amount of time vertex is 
occupied by a type 2 particle grows exponentially with the size of the patch. 

Lemma 3.5 There exists p\ > that only depends on L such that 

P (nf(0) = 2 for all t£(nT + 1, nT + 2) | (0, n) is a good site) > 2p x . 

Proof. Let C be the (finite) set of the configurations restricted to So such that (0, n) is a good site 
and, for each configuration rj £ C, let E v be the set of the realizations of the graphical representation 
restricted to the space-time box Bq x (nT, nT + 2) such that 

{VnT = V on B } n E v C {r?f(0) = 2 for all t G (nT + 1, nT + 2)}. 

Since the events {fj^j, = r\ on So} and T 1 ^ are independent, we have 

P (r?f (0) = 2 for all i G (nT + 1, nT + 2) | r^ T = 77 on B ) 

> P(E V \ r]% T = von B ) = P (E n ) > 0. 

Using that C is finite, we obtain 

P (^(0) = 2 for all t G (nT + 1, nT + 2) | (0, n) is a good site) > inf P N {E V ) > 0. 

)7SC 

Finally, since the events T^ are measurable with respect to the graphical representation restricted 
to the space-time box Bq x (nT, nT + 2), the bound inf^c* P (E v ) only depends on L. □ 

Lemma 3.6 Let P2 > Pi > P c and Ik = exp(cK) as in Lemma \3.4[ Also, denote by £ the Lebesgue 
measure on the real line. Then, for any s G (0,Tk) and K large, 

P (£{* G (s,s +1 K ) : r?f (0) = 2} < K^Zk \ t n > 31 K ) < C 18 exp(- 7l8 K) 

for suitable constants C\% < 00 and 718 > 0. 

Proof. The condition tjv > 31k together with Proposition 13.21 implies the existence of an in- 
all-direction expanding region centered at which contains a positive density of good sites. In 
particular, there exists P2 G (0, 1) such that 

Pq (card {n : nT G (s, s + Ik — 2) and (0, n) 

is a good site} < p%T~ Ik I tn > 3X^) < C%% exp(— 719X). 
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From the previous estimate and Lemma 13.51 it follows that, for K sufficiently large, 

P (£ {t £(s,s + X K ) : n?(p) = 2} < K-H K \ t n > 31 K ) 

< P (£{t G (s, a+I K ): 7?f (0) = 2} < Vx p 2 T- x X K \ t n > 31 K ) 

< P (card {n : nT G (s, s + Z K - 2) and 7?f(0) = 2 for all ts(nT+ 1, nT + 2)} 
< PiP2T' 1 1k I Tzv > 3Xx) < Cis exp(-7i 8 K) 

for suitable C\% < oo and 7i§ > 0. This completes the proof. □ 



4. Proof of Theorem [2] 

In the previous section, we proved that the two-scale multitype contact process restricted to a 
single patch viewed at the L-cube level dominates oriented percolation on Qk- In this section, we 
rely on consequences of this result, namely Lemmas 13.41 and 13.61 to prove that the process on the 
entire lattice viewed at the patch (or iV-cube) level dominates, in a sense to be specified, oriented 
percolation on Q. This will prove in particular Theorem [21 First of all, we consider the interacting 
particle system whose state at time t is a function fjt : % d — > {0,1,2}, and whose evolution at 
vertex x £ A z = Nz + Aq is described by the transition rates 

co-+i(x,fj) = fa Y, !{*) = !} + 2dBi 

x~y 

c 1 ^Q{x,ff) = c 2 -> {x,fj) = 1. 

The dynamics are the same as for the process rjt except at the center of the patches in which first 
type 1 particles appear spontaneously at rate 2dBi and second births of type 2 particles originated 
from adjacent patches are only allowed if the patch is void of 2's. Note that 2dB\ is the rate of 
spontaneous births of type 1 in the process introduced in Section [3] but also an upper bound of the 
rate at which the center of a patch becomes occupied by a 1 originated from an adjacent patch in 
the two-scale multitype contact process. In particular, starting from the same initial configuration, 
the processes rjt and fjt can be coupled in such a way that 

{x £Z d : fj t (x) = 2} C {x G Z d : rj t (x) = 2} 

so it suffices to prove Theorem [2] for the process fjt. The process fjt viewed at the patch level will 
be coupled with the oriented percolation process on Q via the following definition. 

Definition 4.1 Let Xk as in Lemma \3.4l Then, site (z,n) G Q is said to be type 2 stable if 

£{t£nl K + (0,l K ):fjt(Nz) = 2} > K~ X T K . 

The set of sites which are type 2 stable at level n is denoted by 



X n = {z G 7L d : (z, n) G Q and (z, n) is type 2 stable}. 
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The following proposition can be seen as the analog of Proposition 13.21 While Proposition 13.21 is 
concerned with the two-scale multitype contact process restricted to a single patch viewed at the 
intermediate mesoscopic scale, Proposition 14.21 is concerned with the unrestricted process viewed 
at the upper scale. Theorem [5] is a straightforward consequence of Proposition 14. 21 

Proposition 4.2 Let 02 > Pi > /3 C - Then, for K = K{e) sufficiently large, the processes can be 
constructed on the same probability space in such a way that 

P (W n C X n for all n > | W = X ) = 1. 

Proof. Since the evolution rules of the process are invariant by translation of vector u G N7j d , it 
suffices to prove that 

P ((ei, 1) is type 2 stable | (0, 0) is type 2 stable) > 1 - e 

for all K sufficiently large. We assume that site (0, 0) is type 2 stable and let 

<r Cl = inf i t > : Yl ifeH 2} = for all s £ t + (0, 31 K ) > . 

^ weA ei ' 

In words, cr ei is the first time a successful invasion occurs, where successful invasion means that 
a 2 originated from an adjacent patch is sent to A ei and its family survives at least 31k units of 
time in the patch A ei . The aim is to prove that P (<r ei > Zr-) is small for K large. This, together 
with Lemma 13.61 will imply that site (ej, 1) is type 2 stable with probability arbitrarily close to 1 
for large enough K. To estimate the random time a ei we let so = and define by induction 

n = inf It > s.-i : l{fjt(w) + 2} = o| 

and Si = w£ lt>n: Hv t (w) ^ 2} = 1 1. 

^ weA ei ' 

In words, rj is the ith time a type 2 originated from an adjacent patch is born at the center of 
patch A ei and Sj the iih time patch A ei becomes void of 2's. By letting 

M = inf {i > 1 : s, - n > 31 K }, 

we obtain <j ei = vm- Let fj = Sj — T{, and let Sj be the amount of time vertex is occupied by a 
type 2 between time and time r^. Since vertex is occupied by a 2 at least K~ 1 Xk units of 
time until Xk (recall that (0,0) is type 2 stable), on the event {r m > Ik}, we have 

m 

^2 ri + Si > K~ X 1 K , 
i=l 

which implies that 

m 

{r m >I K } C (J {max(fj,Sj) > Z K /2mK}. 
i=i 
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Putting things together, we obtain 



P {a ei >Z K ) = J2 P ( r ™ > Z k and M = m) 



m=l 

oo m 



LXJ III / _ 

E E P (^k < n±K,<*n<\M-n,) (22) 



m=l i=l 

oo m 



m=l i=l 



+ E E P (^>i^™dM = m 



We estimate the right-hand side of (f22l) in three steps (see (j23"j) - (f25"j) below). First of all, observing 
that from time r% to time Sj the process 77^ restricted to A ei evolves according to the transition 
rates of the process 77^ and applying the Markov property, we have 

P(M = m\ M>m) = P (M = 1) = P (t n > 31 K ) 

for all 77t > 2, which implies that 

P(M = m) = (l-psr^Ps with p 3 = P (t n >31 k ) > 0. 

In particular, there is a large m e , fixed from now on, such that 

oo m oo 

Y, P ( M = m ) ^ E rnil-pzT- 1 < i (23) 

m=m e i=l m=m e 

In other respects, Lemma 13,41 implies that 

m t m 



m=l i=l v / 



2m f K 



s™^(^<™^ (24) 



C12 exp(-7i 2 i ; sr) + C13 exp 



713^ \ 
2m e i<C J 



< - 

~ 5 



for K sufficiently large. Finally, letting X ~ Exponential (-B2) and using that Ne± is in state a 
fraction of time of less than K~ l with probability less than e/5m 2 for K large, we obtain 



m=l i=l 



y y P 5 . > < m 2 P 5 > 



^ 6 , 2p(y v ^ U f 2 / \ , 2e 

< - + ttt P \ X > rr < - + mz exp ^ < — 

~ 5 V 2m t K 2 J ~ 5 e V 2m e ET 2 / ~ 5 

for if large. Applying Lemma 13.61 with s = Ik — cr ei , we also have 

P (£ {t £ (1 K , 21 k ) ■■ m(Nei) = 2} < K~ X T K \ a ei < I K ) 

= P ((ei, 1) is not type 2 stable | a ei < Ik) < - 



(25) 



(26) 
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for K sufficiently large. From (123j) -(|26j). we conclude that 

P((ei, 1) is not type 2 stable | (0,0) is type 2 stable) 

< P (a ei > Ik) + P ((ei, 1) is not type 2 stable and a ei < Ik) < e 
for K sufficiently large. This completes the proof. □ 

5. Proof of Theorem [3] 

In this section, we prove that the parameter region in which type 1 and type 2 coexist for the two- 
scale multitype contact process has a positive Lebesgue measure, which contrasts with Neuhauser's 
conjecture about the multitype contact process on the regular lattice. The strategy of our proof is 
as follows. First of all, we establish coexistence for a particular point of the space of the parameters 
by comparing the process with 1-dependent oriented percolation. Interestingly, survival of type 1 
and survival of type 2 are proved by considering different time scales. In other words, the process 
will be simultaneously coupled with two different oriented percolation processes, one following the 
evolution of type 1 particles at a certain time scale, the other one following the evolution of type 2 
particles at a slower time scale. The suitable time scale for type 2 is fixed afterward and depends 
on the time scale chosen for type 1. In both cases, however, the process is viewed at the same 
spatial scale, namely the upper mesoscopic scale (patch level). Since our proof relies on a block 
construction, standard perturbation arguments imply that the coexistence region can be extended 
to an open set containing the coexistence point, which proves Theorem [3l In order to compare the 
process with oriented percolation, we introduce the following definition. 

Definition 5.1 Let Xk as in Lemma \3.4\ and J > 0. Then, site (z,n) G Q is said to be 

1. type 1 stable whenever vertex Nz is occupied by a 1 at time nj . 

2. type 2 stable whenever £{t G nl K + {0,1 K ) : r) t (Nz) = 2} > K- l X K - 
For i = 1,2, the set of type i stable sites at level n is denoted by 

X l n = {z £ Z d : (z,n) G Q and (z,n) is type i stable}. 

Note that the definition of type 2 stable is slightly different from the one in Definition 14.11 in that 
it now applies to events related to the two-scale multitype contact process instead of the modified 
process introduced in Section However, Proposition 14.21 still holds for X% since the set of 2's in 
the two-scale multitype contact process dominates the set of 2's in the modified process. To exhibit 
a point of the space of the parameters at which coexistence occurs, we fix 

Pi > P 2 > ft = fa > /3 C 62 = 1. (27) 

The condition 82 = 1 is to fix the time scale. The condition f3% = indicates that type 1 particles 
can only survive by jumping from patch to patch. In particular, to prove that they survive, the idea 
is to choose 5\ > so small that a 1 at the center of a patch can produce and send its offspring to 
adjacent patches before being killed. More importantly, since f}\ = 0, survival of type 1 particles 
does not depend on the patch size. Coexistence is then obtained by choosing N so large that type 2 
particles can establish themselves an arbitrarily long time in a single patch. Since type 1 particles 
have a positive death rate, centers of patches are empty a positive fraction of time which allows 
type 2 particles to survive by invading adjacent patches from time to time. 
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Proposition 5.2 Assume <\27\i . Then, for suitable J < oo, 5\ > and N < oo, the processes can 
be constructed on the same probability space in such a way that 

P (VP* C Xl and W% C X 2 n for alln>0\ = and W$ = X$) = 1 

where and W% are two copies ofW n . 

Proof. As in Proposition 14.21 it suffices to prove that, for i = 1,2, 

P ((ex, 1) is type i stable | (0, 0) is type i stable) > 1 — e (28) 

for a suitable choice of the parameters. The first step is to show that there is enough room for type 
1 particles to invade the center of the patch A ei . By observing that 

c ^ 2 (Ne 1 ,7 1 ) < 2d{B 2 + (3 2 ) and 02^0(^1,77) = S 2 = 1 

we obtain 

hmsup E[s- 1 2{te(0,s): Vt (Ne 1 ) = 2}] < G 2 := ; < 1- 

s-i-oo 1 + la [B 2 + p 2 ) 

Let eo 6 (0, 1 — 62). Large deviation estimates for the Poisson distribution give 

P (£ {t G (0, J) : Vt(N ei ) = 2} > (0 2 + e ) J) < C 20 exp(- 720l 7) 
for suitable C20 < 00 and 720 > 0. In particular, 

P (rjtiNex) ^ 1 for all t G (0, J) and »ft(0) = 1 for all t G (0, J)) 

< P (£ {t G (0, J) : %(iVei) = 2} > (0 2 + e )J) 

+ P (r? t (iVei) / 1 for all t G (0, J) and r/ t (0) = 1 for all t G (0, J) and 
£{te (0, J) : Vt(Nex) = 2} < (0 2 + e )J) 

< C 20 exp(- 7 2oJ) + exp(-B 1 (l-Q 2 -e )J). 

Finally, taking J large and then S\ > small, we get 

P((ei, 1) is not type 1 stable | (0,0) is type 1 stable) = P{r ) j(Ne l ) / 1 | 770(0) 

< P (r]j(Nei) / 1 and r) t (Nei) = 1 for some t G (0, J7") | 770(0) = 1) 

+ P (r?t(iVei) / 1 for all t G (0, J) and 7? t (0) = 1 for all t G (0, J) | 770(0) 
+ P(f]t(0) / 1 for some t G (0,J) | 770(0) = 1) 

< C 20 exp(-72 0l 7) + exp(- J B 1 (l-e 2 -e ) l 7) + 2 (1 - exp(-5 ll 7)) < e 

which establishes (|28p for i = 1. Moreover, survival of type 1 particles holds regardless of the 
patch size so the proof that conditions (j28|) for i = 1 and 7 = 2 hold simultaneously for the same 
parameters follows by taking L and if sufficiently large, and applying the results of the previous 
two sections. This proves that both types coexist. To conclude, we briefly justify the fact that the 
results of the previous sections hold as well under the new assumptions 



Pi = and Si > 0. 



(29) 
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First, the condition /3i = implies that, starting from any initial configuration, 

p (rjf(x) = 1 for some x 6 B z \ {0} and t > VI) < L d exp(-£iVT) 

which, now that Si > is fixed, can be made arbitrarily small by taking L large. In words, except 
at the center of the patch, all the l's in B z C Aq are rapidly killed, so the proof of Lemma 13.31 
extends easily under (129f) . Proposition 13.21 and Lemma 13.41 follow as well. Now, we observe that 
the condition 5i > implies that vertex is empty a positive fraction of time, which allows for 
invasions of particles of type 2 at the center of the patch. With this in mind, one can check easily 
that the proofs of Lemmas 13.51 and 13.61 also apply under the assumptions (j29[) . Note however that 
the new lower bound p\ > in Lemma 13.51 might be smaller. Finally, the proof of Proposition 14.21 
still holds as a consequence of Lemmas 13.41 and 13.61 □ 
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